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Abstract 



^vq An important tool in the field of metamaterials is the extraction of effective material parameters from simulated or measured 
y—{ scattering parameters of a sample. Here we discuss a retrieval method for thin-film structures that can be approximated by a two- 
<^> dimensional scattering sheet. We determine the effective sheet conductivity from the scattering parameters and we point out the 
importance of the magnetic sheet current to avoid an overdetermined inversion problem. Subsequently, we present two applications 
O ,of the sheet retrieval method. First, we determine the effective sheet conductivity of thin silver films and we compare the resulting 
D conductivities with the sheet conductivity of graphene. Second, we apply the method to a cut-wire metamaterial with an electric 
dipole resonance. The method is valid for thin-film structures such as two-dimensional metamaterials and frequency-selective 
\Q surfaces and can be easily generalized for anisotropic or chiral media. 

(N 
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1. Introduction 

The development of metamaterials, i.e., artificial microstruc- 
tured materials in which small, subwavelength electric circuits 
replace atoms as the basic unit of interaction with electromag- 
netic radiation, has been a vibrant research topic in the field 
of nanophotonics during the past decade Q3E1E1. Structur- 
ing metamaterials on a subwavelength scale makes it possible 
to create electromagnetic media with properties not found in 
natural materials, while still allowing to describe them as ef- 
fectively continuous media with constitutive parameters such 
as the electric permittivity and the magnetic permeability PI . 
With appropriately designed constituents, it was shown feasi- 
ble to design metamaterials with exotic material response, e.g., 
magnetism at terahertz and optical frequencies, simultaneous 
negative permittivity and negative permeability (the so-called 
left-handed materials) 0, giant chirality [6|, and slow-light 
media [7 8, 9|. The first metamaterial was fabricated in 2001 
by combining metal wires exhibiting negative permittivity and 
split-ring resonators (SRR) exhibiting negative permeability in 
a single material HOI . SRRs are still commonly used in the 
microwave band, but have been replaced by other magnetically 
resonant structures such as slab-wire pairs and fishnets at tera- 
hertz frequencies and above, since these meta-atoms ease the 
problem of saturation of the magnetic response at those fre- 
quencies Ifm and also simplify the fabrication. 

Metamaterials exhibit novel electromagnetic phenomena, 
such as backward wave propagation, negative refraction, and 
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inverse Doppler effect Ifl2l . Furthermore, they enable the com- 
pensation of the propagation phase and the restoration of the 
amplitude of evanescent waves in optical structures, result- 
ing in lenses with subwavelength resolution lfl3ll and photonic 
devices going beyond the diffraction limit iTBl [T5l fl6l H71. 
Through the technique of transformation optics, metamaterials 
with arbitrary values of the permittivity and permeability ulti- 
mately allow for maximal control over light propagation, which 
has led to surprising new physics such as invisibility cloaks, 
beam transformation, and frequency conversion in linear me- 
dia d [H |2Q] HQ . 

A basic tool in the study of metamaterials is the so-called re- 
trieval method, i.e., the extraction of effective medium param- 
eters corresponding to a metamaterial with given microscopic 
structure. Effective material parameters are important because 
they make the link between the microscopic response of meta- 
materials and the homogeneous macroscopic media assumed in 
many proposed applications. We assume here that the meta- 
atoms are sufficiently subwavelength, so that the effective re- 
sponse is only weakly spatially dispersive and, hence, can be 
described by frequency-dependent polarization and magnetiza- 
tion fields [ 22 1 . For ease of presentation, we assume here (two- 
dimensional) isotropy in the plane of the thin-film structure, but 
the considerations below can be straightforwardly generalized 
for anisotropic samples. 

There are several methods to determine the effective material 
parameters. One approach involves the averaging of the elec- 
tromagnetic fields following their definition [23 1 ; the volume 
averaging can sometimes be replaced by line and surface aver- 
ages [ 24- 1 . These averaging methods and other methods requir- 
ing the knowledge of the microscopic fields [25 1 are, however, 
often difficult to apply, since it is complicated or even impossi- 



Preprint submitted to Elsevier 



September 27, 2012 



E inc exp[i(k// • r + kj_- r - cot)] 



RE mc exp[i(k// r - k ± r - cot)] 



a Je 



rE inc exp[i(k//- r + k x - r - cot)] 



Figure 1 : Scattering of a plane wave by a two-dimensional current sheet with 
electric sheet current j e and magnetic sheet current j m . 



ble to obtain the microscopic fields. Another approach involves 
comparing the metamaterial sample with a slab of a homoge- 
neous medium and determining the constitutive parameters of 
the slab that gives the same scattering matrix as the metama- 
terial E6l l27l l28l . This is attractive, because scattering exper- 
iments can be easily set up, both experimentally and numeri- 
cally. One problem that arises when using the latter retrieval 
method with thin-film samples is choosing a suitable length of 
the slab 11291 . Since the thickness of such thin-film samples is 
typically much smaller than the wavelength, it makes actually 
more sense to compare it to a two-dimensional current sheet and 
characterize the sample by sheet susceptibilities |29l or sheet 
currents (this paper) that result in the same scattering parame- 
ters as obtained for the metamaterial. 



2. Retrieval method 

The two-dimensional retrieval method starts from a rather 
standard problem of electrodynamics — the scattered fields of 
a two-dimensional sheet carrying electric and magnetic sheet 
currents. The problem is sketched in Fig. [T] A plane wave 
E; nc exp[i(k a ■ r + kj_ • r - cot)] illuminates the current sheet and 
is scattered into a reflected wave /?E; nc exp[i(k// -r-k ± -r-cot)] 
and a transmitted wave T E mc exp[i(k// ■ r + k ± ■ r - cot)]. Since 
the current sheet is assumed to have a linear response, we 
can work in the harmonic regime, and we will drop the factor 
exp(-iwf) from now on. 

The incident and scattered waves must satisfy the boundary 
conditions 



n (D-D')=p e , 
n (B-B')=p m , 
n x (E - E') = -j n 
nx(H-H')=j e , 



(1) 
(2) 
(3) 
(4) 



where n is the surface normal of the current sheet, j e and j m are 
the electric and magnetic sheet currents, respectively, D is the 
electric displacement field, B is the magnetic induction field, E 
is the electric field, and H is the magnetic field. It might seem 
strange that we include a magnetic sheet current here, but we 



will see that this is essential to the retrieval method and we will 
explain below how an effective magnetic sheet current origi- 
nates in a thin-film sample. 

Substituting the incident and scattered fields in Eqs. Q-Q, 
and using the identity H = k x E/(wjUo), yields 



j e = — (1 - fl-r)nx(kxE inc ) 
cono 

= r 1 a-/?-7 , )(E iIlc ) // , 

j m = -(l+ft-r)nxE mc 

= ^'(l+ J r?-7')(H inc ) // , 



(5) 
(6) 



where £ is the wave impedance defined by = Ic ± /(colio), 
rjo is the characteristic impedance of vacuum, and Icq = co/c is 
the free-space wavenumber. 

To find the conductivities, we have to relate the sheet currents 



to the local fields, j e = cr^ Ei oc and j m = o^""Hi 0C . Since the 
current sheet carries electric as well as magnetic sheet currents, 
both the electric and magnetic fields have discontinuities at the 
current sheet. Therefore, the local fields must be defined as an 
average across the current sheet: 

E/i + W i 
Eioc = = - (1 + R + T) (Eh,,.)// , (7) 

H//+H' i 
H loc = _J/ = _(!_* + DfHteV (8) 

Combining Eqs. |5])-(|6| with Eqs. (j7]i-([8]), we arrive at the for- 
mulae that give us the effective conductivities that are com- 
mensurate with given (measured or simulated) scattering co- 
efficients: 



(e)= 2 1-R-T 
°~U £ \ 1 + R + T 



II h \l-R + T 



l+R-T 



(9) 
(10) 



This is the central result of this paper. 

It is interesting to note that these formulae can be uniquely 
inversed, 



T = 



(e) (m) 
r (T 

II II 



4 + 2^,7 +2£" V/ 



2(^>-rM?) 



a . <~*r (e) . (m) , (e) (m) ' 

4 + 2£o-)[ + 2£ l cr)/ + cr) / 'cr) / 



(11) 
(12) 



ir 11 



from which it is obvious that there is an isomorphism between 
the scattering parameters and the sheet conductivities. Conse- 
quently, the thin sheet retrieval does not suffer from the am- 
biguities encountered in the traditional retrieval method, such 
as branch cuts, residue classes, or free parameters (the slab 
thickness). We can now also understand why we need the mag- 
netic sheet current; without the magnetic sheet current, the re- 
trieval problem would in general be overdetermined with only 
one complex parameter (the electric sheet conductivity) deter- 
mined by two complex scattering coefficients (R and T). Of 
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Figure 2: (a) Effective electric sheet conductivity of thin silver films as a func- 
tion of film thickness at A = 1.55 |am. (b) Effective magnetic sheet conductivity 
of thin silver films as a function of film thickness at A = 1 .55 pm. 



course, if the sample can really be considered as an electric 
current sheet, the reflection and transmission coefficients are 
related by 1 + R = T and the inversion would not be overde- 
termined. Nevertheless, it may probably be better in practice 
to use Eqs. |9]l-( 10 1 and then check whether the magnetic sheet 
conductance is negligible. 



Figure 3: Unit cell of the single-layer cut-wire metamaterial considered here; 
the dashed line indicates the unit cell boundaries. The wire has large end caps 
to keep the lattice constant sufficiently below the wavelength. 



resulting sheet conductivities are plotted in Fig. |2ja)-(b). The 
magnetic sheet conductivity is negligible for the thinnest films, 
but its imaginary part reaches 1 % of the vacuum impedance at 
a thickness of approximately 25 nm. For thicker films, the re- 
sponse of the silver films can no longer be described properly 
by a sheet model. A 25 nm silver film has an electric sheet con- 
ductivity of cry, = 0.48/?7o = 1262 \iS, to be compared with the 
sheet conductivity of graphene equal to cr g = 60.8 |iS. So even 
for relatively thin films, the current transport in silver films is 
more efficient than in graphene at infrared and optical frequen- 
cies. 



3. Application to thin silver films and comparison with 
graphene 

Thin conducting films are frequently employed in nanopho- 
tonics. Gold and silver are normally used because of their mod- 
erate losses and relatively easy chemical handling. Neverthe- 
less, dissipative losses in those metals are still appreciable at 
optical frequencies and it is, therefore, interesting to search for 
alternative conducting materials. Here we want to compare a 
silver film with graphene. The conducting properties of silver 
are characterized by its frequency-dependent bulk conductivity 
in S/m. As a two-dimensional electronic system, graphene is 
characterized by its sheet conductivity, which is almost flat and 
equal to <x g = 6.08 x 10~ 5 S starting in the infrared and extend- 
ing into the visible frequency band. 

In order to compare the current transport in thin silver films 
and graphene, we use our sheet retrieval method to describe 
the properties of silver films with thickness between 10 nm and 
40 nm by sheet conductivities. We first calculate the scattering 
parameters for silver films at a wavelength of 1.55 [J.m using 
the Drude model from Ordal et al. l30ll for the bulk conduc- 
tivity (note that the Drude model in Ordal et al. has been fit- 
ted against experimental data from thin films consistent with 
the film thickness here) and then we apply Eqs. (|9]l-([T0]>. The 



4. Application to a single-layer metamaterial 

As a second example, we determine the effective sheet con- 
ductivity of a single-layer metamaterial. The unit cell of the 
metamaterial, shown in Fig. [5] consists of a simple cut-wire 
made from copper. The cut-wire has a quasistatic electric 
dipole resonance. We have calculated the scattering parame- 
ters using a frequency-domain electromagnetics solver (CST 
Microwave Studio). One unit cell is simulated with mirror sym- 
metry boundary conditions (PEC/PMC) to model a periodic ar- 
ray of cut-wires. The resulting scattering parameters and the 
absorption of the metamaterial are plotted in Fig. [4] We see that 
there is a resonance — the electric dipole resonance — at 25 GHz 
where the incident wave is highly scattered and we observe a 
peak in the absorption around the same frequency. 

The retrieved electric sheet conductivity following Eqs. (|9|l- 
(T0) is plotted in Fig. [5] (The magnetic sheet conductiv- 
ity is several order of magnitudes smaller than the electric 
sheet conductivity.) The electric sheet conductivity has a clear 
Lorentzian resonance, a fact that we have checked by fitting a 
Lorentzian function, 

oM = —™l (13) 
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Figure 4: Scattering parameters and absorption of the single-layer cut-wire 
metamaterial considered here. A resonance around 24 GHz is observed. 

to the complex conductivity (the modification of the Lorentzian 
is because the current is the time derivative of the dielec- 
tric polarization). The Lorentzian spectrum fits the retrieved 
conductivity almost perfectly with resonance frequency / r = 
24.3 GHz, damping frequency F = 0.36 GHz, coupling con- 
stant k = 88.3 MHz, and background polarization p - 1.89 fs. 
The excellent fit of this Lorentzian spectrum to the electric con- 
ductivity confirms the quasistatic nature of the electric dipole 
resonance. The latter fact, together with the negligible magnetic 
conductivity, confirms that a single-layer cut-wire metamaterial 
can be almost perfectly described by the sheet model. 

5. Conclusions 

In this paper, we have detailed a retrieval approach for thin- 
film systems by mapping the scattering parameters (reflection 
and transmission) onto the electric and magnetic sheet conduc- 
tivities of a two-dimensional current sheet. We have applied the 
sheet retrieval method to two nanophotonic systems. For thin 
films, we find the sheet retrieval to be valid for films with thick- 
nesses up to 25 nm at A — 1.55 u.m and we find the sheet con- 
ductivity of these silver films to be larger than the conductivity 
of graphene. For the single-layer cut-wire metamaterial, we 
demonstrate that the effective sheet conductance is an excellent 
representation of the scattering response. We believe that the 
sheet retrieval may be more appropriate for many single-layer 
metamaterial s, because it does not require any model parameter 
(e.g., the thickness) to be determined. 
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